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Abstract 

We study the decay of “false” domain walls, which are metastable states of the quantum theory 
where the true vacuum is trapped inside the wall, with the false vacuum outside. We consider a 
theory with two scalar fields, a shepherd field and a field of sheep. The shepherd field serves to 
herd the solitons of the sheep field so that they are nicely bunched together. However, quantum 
tunnelling of the shepherd field releases the sheep to spread out uncontrollably. We show how to 
calculate the tunnelling amplitude for such a disintegration. 

PACS numbers: 11.27.+d, 98.80.Cq, ll.15.Ex, ll.15.Kc 


address after Sept. 1st, 2015: School of Physics and Astronomy, University of Minnesota, 116 Church St 

SE, Minneapolis, MN 55455 USA 
* M.Haberichter@kent.ac.uk 
1 richard. mackenzie@umontreal .ca 
paranj@lps.umontreal.ca 
^ klingon_ecology@hotniail.com 


1 



I. INTRODUCTION 


The notion of a vacuum in field theory is quite a bit richer than that in everyday language. 
In a theory with scalar fields, any state that is a local minimum of the potential energy 
density is a vacuum. More than one vacuum can exist, and not all vacua need have the 
same energy. The lowest such state is called the true vacuum; any higher-energy local 
minimum is a false vacuum. 

The most common reason for a multiplicity of vacua is spontaneous symmetry breaking. 
Breaking of a discrete symmetry gives rise to discrete vacua whereas breaking of a continuous 
symmetry gives rise to a continuum of vacua. In this paper we restrict ourselves to discrete 
vacua. Depending on the model, a number (possibly infinite) of vacua of any energy can 
exist. Field configurations that interpolate between different degenerate vacua, true or false, 
give rise to domain walls (see for example [1]), which themselves are stable or metastable. 
Domain walls are very important in cosmology [2], and in condensed matter physics [3], 
where they quite often arise in the study of first order phase transitions. A full review 
of domain walls and their applications is beyond the scope of this article, and we refer the 
reader to the cited literature. It is clear that domain walls are very important in wide swaths 
of theoretical physics. 

In this paper we consider the special case of domain walls interpolating between distinct 
false vacua, with true vacuum in the core of the domain wall. The very existence of domain 
walls in our model is due to a subtle interplay between the energetics of two fields. The 
domain wall we consider entraps several quanta of a quantum field, 0, which we call a 
field of sheep, which on their own would separate to infinity, leaving behind regions of true 
vacuum. It is the sheep field that is in the false vacuum outside the domain wall, but in 
regions passing through the true vacuum inside the domain wall. These sheep are herded 
into staying together, inside the domain wall, by a field we call the shepherd field, "0, which is 
in its true vacuum outside the domain wall, but in fact in its false vacuum inside the domain 
wall. The shepherd field is unstable to quantum tunnelling to its true vacuum; once this 
occurs, the sheep are without a shepherd and will spread out to infinity. We compute the 
amplitude for such a decay in a specific model analytically, aided by numerical calculations, 
within a well-defined approximation scheme. We have also studied an analogous model [4], 
but in that model, there are not multiple quanta of one field trapped inside the domain wall. 
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II. THE MODEL 


Consider the model defined by the Lagrangian density 


C = - V (^, 0) 


( 1 ) 


with 


V(^, 0) = V4^) + V4</>) + 0) - Co 


( 2 ) 


where 


Ca(V’) = + afid’ - aY + ej)) 

1/^(0) = /3(sin^(7r(^) + e^sin^(7r0/2)) 


( 3 ) 


( 4 ) 


while 


0) = A 


{ip - afUij + af + el)) 
(W)-14(1/2)2 + 7^ 


( 5 ) 


A is considered to be very small compared to a and 7 is considered small compared to 
14(1/2). We impose that V{ijj,(j)) > 0. The first three terms are positive semi-definite, and 
then Vq is chosen so that V{ip, 4>) vanishes at its global minimnm. The ep field is essentially 
a sine-Gordon field (with alternating trne and false vacna) while the ip field has a rather 
standard donble well potential with a small asymmetry. Both fields admit kink type solitons 
when 60 = = 0. The interaction between the two fields is constrncted to disallow the 

passage of the (p solitons throngh the ip solitons. The specific potential that we use is not 
important; the only feature that is crucial is this “non-commutativity” of solitons, and any 
potential that achieves the result is adequate. 
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FIG. 1. (Color online) The potential for the field ip with a = 1. 
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FIG. 2. (Color online) The potential for the field cj). 

It is easy to see that the global minimum of occurs at V’ = —cl while a local minimum 
occurs near -0 = a, as can be seen in Fig. (1). These values are mildly perturbed by 4>) 

but for small values of A; however the changes can be made arbitrarily small. The minimum 
of occurs at even integers (j) = 2k for any integer fc, with a slightly higher local 

minimum occurring at odd integers (p = 2k + 1, and local maxima occur at each half odd 
integer, as can be seen in Fig. (2). Again, these values are mildly perturbed by (p) but 

for small values of A the changes can be made arbitrarily small. 

Obviously for A = 0, the true vacuum corresponds to {'ip = —a, (p = 2k} while {pj = 
—a, (p = 2k + 1} or {pj = a, (p = 2k or 2k + 1} correspond to the false vacuum. Of course 
with A non-zero, the values of the helds at the minima will be slightly modihed, however 
the structure of the false and true vacua will remain unchanged. In the ensuing discussion 
we will identify the vacua by the values that the helds take when A = 0, but always with 
the understanding that the actual values that the helds take are slightly perturbed. 

As there are several discrete vacua, we can easily imagine conhgurations corresponding to 
domain walls. Consider a domain wall interpolating between the two false vacua, (p = 2k+ 1 
and (p = 2k' -\- 1 but pj = —a, its true vacuum. Thus (p interpolates from 2/c -|- 1 to 
2k' + 1 while the held ip = —a remains at its true minimum. Such a domain wall will 
be immediately unstable to separating into, in principle, several subsidiary domain walls. 
Consider the specihc example where the (p held interpolates from 0 = —1—)-0 = l. Such 
a domain wall will be unstable, separating into two half domain walls that interpolate from 
0 = —1—>-0 = 0 and then from 0 = 0 —)■ 0 = 1 while the ip held remains always at 0 = — a. 
The intermediate region being true vacuum, 0 = 0,0 = —a, the two half domain walls 
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feel a repulsive force and will accelerate away from each other leaving behind them true 
vacuum. However, taking into account the interaction term, such a false domain wall can be 
rendered metastable (classically stable, but unstable to quantum tunnelling) by imposing 
some transitions in the 'ip held to its false vacuum. The pj held must makes its transition 
to the left of the point where (p held begins its transitions, and ip must return back to its 
true vacuum only after (p has completed its transitions and has reached 1. There would be a 
soliton and an anti-soliton pair in the ip held. The conhguration corresponds to something 
like Fig. (3), which is a numerical solution of the equations of motion, and has been obtained 
using the gradient how method [5]. 

In the stability analysis, the interaction term makes the following crucial contribution. 
The point is that if (p passes through the value 1/2 when ip is such that the numerator in 
the interaction term is not very small or even zero (this means that ip ^ ±a), then the 
denominator of the interaction term becomes 7 ^, which is chosen to be very small, and thus 
the interaction term can be made to blow up. In this way, the contribution to the energy 
of the interaction term can be made arbitrarily large, which we assume, and therefore such 
a situation will be energetically disfavoured. What this means specihcally is that the (p 
solitons may not pass through the ip solitons (in which necessarily occur regions where the 
ip held is signihcantly different for ip = ±a) without the energy blowing up. There is an 
energy barrier to such a process. Thus if ip is forced to make a transition from —a to a only 
before or after cp makes its transitions, the solitons of the (p held cannot pass through those 
of the ip held without encountering an insurmountable energy barrier. But then the solitons 
of the (p held, the sheep, are trapped within the solitons of the ip held, the shepherd. 

Inside the domain wall, two adjacent (p half-solitons (say, one going from (p = 2k to 2k + l 
and the other from 2k + 1 to 2k + 2) attract each other. This is because the energy is 
minimized by reducing the region over which <p is in its false vacuum, (p = 2k + 1. Thus the 
net (p soliton passes from even integer to even integer. But then these full solitons further 
minimize their energy by separating as far as possible from each other, as this minimizes any 
gradient energies that occur due to their proximity; hence, they have a repulsive interaction. 
The half-solitons on either end are also simply repelled by the adjacent full solitons. 

Thus in the example shown in Figure (4), in which (p interpolates from -3 to 3, the 
half-soliton from -3 to -2 is repelled by the full solitons from -2 to 0, which also repels the 
subsequent full soliton from 0 to 2 which in turn repels the hnal half-soliton from 2 to 3. 


5 


The il) (shepherd) solitons have the false vacuum 'ip ^ a in between them, so they attract one 
another. The sheep want to separate as far as possible from one another, but they cannot 
do so without pushing the shepherd solitons further and further apart. The energy in doing 
so, increases linearly with the separation of the shepherd solitons. The competing energy 
decrease coming from the separating sheep is at best decreased to a constant value (the total 
mass of sheep) as the separation between the sheep becomes large. Hence, energetically, the 
shepherd solitons force the sheep solitons to be herded together at a finite size where the 
outward pressure from the sheep is balanced against the inward vacuum pressure pushing 
the shepherds together. The whole system forms a classically stable configuration. An even 
more extreme example is given in Figure (5) and the corresponding zooms in Figures (6), 
which contains 31 sheep. 



FIG. 3. (Color online) The metastable domain wall configuration, cp passing from —1 —>■ 1 and ip 
passing from — 1 —)• 1 —)■ — 1. For this and all subsequent figures, the parameters were chosen as 
a = 1, 60 = 0.01, 60 , = 0.5, A = 0.1, a = 0.5, /3 = 0.5 ,7 = 0.01 

However, quantum mechanically, the herded sheep are actually in a metastable config¬ 
uration. Tunnelling processes can permit the pj field to make a quantum transition to its 
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FIG. 4. (Color online) The metastable domain wall configuration, (j) passing from —3 —)• 3 and ^|J 
passing from — 1 —?■ 1 —)• — 1. 

true vacuum at all points where the (j) field is near its true vacuum (and specifically away 
from values where = V{l/2)). Such transitions truly liberate the sheep, which can 

now move apart freely. 

III. ANALYTICAL CALCULATION OF THE DECAY RATE IN THE FREE KINK 
APPROXIMATION 

As we have just noted, the soliton is metastable: it will decay via quantum tunnelling. 
The tunnelling transition is mediated by an instanton, which is called a bounce. The bounce 
corresponds to a trajectory in Euclidean time, t —)■ —ir. This changes the sign of the kinetic 
energy, T —)■ —T, and therefore bounce configuration can be thought of as motion in the 
potential —V where we include the spatial gradient terms as part of the potential. The 
Euclidean action can be written as Se = T-hV and the conserved Euclidean “energy” is E = 
T—V. The bounce begins at the metastable configuration, arrives at the configuration which 
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FIG. 5. (Color online) The metastable domain wall configuration with 31 sheep, (j) passing from 
— 15 —>■ 15 and -0 passing from —1 —)> 1 —>■ —1. 

will materialize upon the event of the quantum transition, and then “bounces back”, that 
is, returns to the initial configuration. Therefore, the bounce will begin at the metastable 
soliton as in Fig. (7a) followed by an evolution during which the field V’ which is initially 
equal to a in the region in between the two half-solitons sinks down until it finally makes it 
to "0 = —a, corresponding to a configuration of two free half-solitons as depicted in Fig. (7b). 
This configuration would actually correspond to the unstable configuration at the top of the 
potential barrier. The two half-solitons must be created, however, with some kinetic energy, 
since during the Euclidean evolution, the “energy” T — 17 is conserved, and —17 is more 
negative at this configuration of two essentially free solitons. In terms of the Minkowski 
dynamics, since this configuration corresponds to the top of the barrier, it has much more 
potential energy, -|-I7, than the initial configuration. Then in the Euclidean dynamics, it has 
potential energy —17; hence there has to be compensating kinetic energy. The solitons then 
separate and rise up the (negative) potential until they stop, after which they bounce back. 
When they reach the bounce point, that is the configuration at which they will materialize 
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FIG. 6. (Color online) (a) Zoom of the metastable domain wall confignration with 31 sheep for the 
held (p. (b) Zoom of the metastable domain wall conhgnration with 31 sheep for the held ip. 


when the quantum tunnelling transition occurs, and from that point on Minkowski evolution 
takes over. 

The Euclidean Lagrangian is given by 
Le = j dx 

The fields depend on time only through their dependence on the positions of the solitons. 
Indeed, in a first approximation, the fields depend independently on the positions of the 
solitons, which is the appropriate collective coordinate describing their evolution. This 
means we can approximate (f){x,xi,X 2 ) = <pi{x — xi) + (j) 2 {x — X 2 ) while 'ip{x,xi,X 2 ) = 
■01 (x — xi) +'ip 2 {x — X 2 ) + 'ipQ where -00 is a necessary constant offset. This approximation is 
certainly valid when the solitons are separated from one another. It is not valid when they 
are very close, and we will have to resort to numerical calculations to find the action there, 
but we leave that for later. For example, for 'ip we get 


^ 0 ^ + ^ 0 ^ + ^ 0 '^ + ^ 0 '^ + V ( 0 , 0 ) 


( 6 ) 


0 = d^,'lp(x,Xi,X2)Xi + d^2'lp(x,Xi,X2)x2 

= d^.-lpiix - Xi)xi + 0a:202(a^ ” X 2 )X 2 

=-d^-ipiix - Xi)xi - d^i’2(x - X2)X2. (7) 


dx'iplix — xi) = 0( is a function that is sharply peaked around xi while dx'ip 2 {x — X 2 ) = 02 
is sharply peaked around X 2 ; thus, the product of these functions essentially vanishes. This 
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gives 



( 8 ) 

Furthermore, the nature of the solitons at xi and X 2 in each held are essentially identical; 
they are simply the “kink”-type solitons of an (almost) doubly degenerate (actually multiply 
degenerate, but only adjacent minima are relevant here) potential. Therefore we will com¬ 
pute the contribution of the kinks in the limit that the potential is exactly degenerate, and 
we note that kink or antikink obviously give an equal contribution. Thus we write 


j dx^cj)'^ = j dx^<p 2 = ^ and y dx^ij'^ 


dx-ip 



( 9 ) 


and should not be confused with the actual masses of the perturbative excitations 
of each respective held. The contribution of the potential to the action will also be easily 
deconstructed. There will be a contribution from around the region of each soliton and there 
will be a contribution from the region in between the solitons when they separate. We hnd 
that the interaction potential plays no role in evaluating the action. We will assume this for 
the moment, and discuss it in more detail in the next subsection. 


A. Free kink approximation 

We can write 


dxV{(t),tp) 


dx iVfj) -T V^) dx (V(j) V^) 


( 10 ) 


' X'^Xl 


' X'^X2 


and we will compute each term in Eqns.(9,10) in the free kink approximation. For a held y, 
which satishes an equation with solution given by a free kink, we have a Lagrangian 

1 


L = -y'^ + E(y) 


and equation of motion 


x" - v\x) = 0. 

The hrst integral of the equation of motion is conservation of energy: 

\x‘ - V(x) = c ^ 0 


( 11 ) 


( 12 ) 


(13) 
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where C is a constant which we normalize to zero. Therefore x' = ■\/21/(x). Then the 
integrals that appear in the action above are rednced to qnadratnre 


rdx^-x'^= r dx^-xmx)x'^ r dxJ^ (w) 

'J Xi Xi j 

and 

£’ dxVix) = £' dx^^2V{x)x' = dx\f^; (15) 

evidently, the integrals are eqnah 


B. Analytical evaluation of the Euclidean action 


For 0, the potential is = Q;sin^(7r0) which gives 


^ = r = 


dxV^{(j)) = \ - d(j){- sin(7r0)) = W - 




a cos{7i(j)) 


71 


\Pl(X 


-1 


TT 


(16) 

The minns sign in front of the sin(7r0) comes from the sqnare root since for the first kink 
(j) G [—1, 0] where sin(7r0) < 0. For the field '0 we have = (3 ("0^ — which gives 
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(17) 


We note the sign change in the intermediate eqnation as ("0^ — n^) is negative for -0 G [—a, a] 
and the limits x, and Xf are assnmed to be both close to Xi or X 2 as reqnired. 

We will evalnate the energy, first in the metastable confignration and second in the 
configuration at the top of the potential barrier, so that eventually we can normalize the 
Euclidean action properly and then evaluate it. For the metastable configuration, the energy 
essentially obtains contributions from the locations of the kink solitons, two for (j) and two 
for ip, hence 


Eo = - dx + p'i + V’f + V'?) + / dxV (0,0) 


XPSX2 


= {m^ + m^)+ dx(l/(0) + E(0)) + / (ix(E(0) + E(0)) 

J X'^X\ J X^X2 

= 2(m0 + m^) 


(18) 


using Eqns. (16,17). We will subtract Eq from the potential so that we normalize the initial 
configuration to vanishing energy. 
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For the configuration at the top of the potential barrier we must take into account of the 
following differences: the solitons must have kinetic energy in order to conserve energy, we 
have two additional kinks in the V’ field, and if the solitons separate by a distance \x 2 — X\\ 
the energy increases (it actually decreases as this contribution is negative) by dxV (0, V’)- 
For the metastable configuration we did not consider that the solitons could separate; it is 
assumed that they are herded together by the shepherd field. If the solitons at the top of 
the potential barrier separate, they leave behind the true vacuum. Since we will normalize 
the vacuum energy so that the false vacuum outside the shepherd soliton has zero energy 
density, the true vacuum between will have negative energy, so this contribution in fact 
decreases the energy. The 'ip field is in its true vacuum outside the shepherd soliton; thus 
there is no energy density difference for it. However, the 0 field is in its false vacuum outside 
the shepherd soliton; hence, looking at the potential Eqn. (4), we see that the contribution 
will be 

rx2 

/ dxV{(l),'ip) =-ae^\x2-xi\. (19) 

J Xi 

Then the energy for the configuration at the top of the barrier will be 

Et - Eo = ^(m^ + 2m^){xl + xj) + 2{m^ + 2m^) - 2{m^ + m^) - ae 4 ,\x 2 - xi| 

= ^(m^ + 2m^){xl + xl) + 2m^ - Q;e</,|x2 - xi| 

= T + V (20) 

The Minkowski action is Sm = j — V) while the energy is E = T + V. Analytically 

continuing to Euclidean time, T —>■ —T and Sm ~ f dx (T -h V) = —Se- Then the 

Euclidean action is Se = J dr {T + V) and is easily read off as 

‘S'e = y dr Q("i</. + 2m^){xl + xl) + 2m^ - ae^\x 2 - (21) 

where the overdot now means the derivative with respect to r. The Euclidean trajectory, 

which we can compute analytically, starts at the configuration at the top of the barrier. 

The solitons separate until they reach the bounce point, and then they come back together 
and coalesce. We define center of mass and relative coordinates X = [xi + X2)/2 and 
X = X 2 — Xl, in terms of these, xl + xl = 2X^ + ^x^. We take the center of mass coordinate 
to be constant, W = W = 0 so that x^ + xl = ^x^. The Euclidean “energy”, T — V, is 
conserved and normalized to zero 

0 = T -V = ^ j - 2m^ + ae^x (22) 


12 


where we restrict a; > 0. The solitons separate until x = 0, hence at that point x = xb 
which is given by 


2m. 


Xb = 




ae^ 


(23) 


Then the contribution to the action from this part of the trajectory will be 

Se = j dr^ x"^ + 2m^ - a6^\x\. (24) 

The energy condition Eqn. (22) implies 


X = 2a 


{2m^ — ae^x) 


(25) 


+ 2m^ 

Replacing, in the usual way as in Eqns. (11-15) for x in the first term of Eqn. (24) we find 

i (? + ’"-'■) ^ + 2m*i 

J dT^2m^ - ae^x^m^ 2m^x 

-2^m^ + 2m.^{2m^ - ae^xY^"^ 


Sk = / dr- 


/ dx^j2m^ — ae^x^m^ + 2m^ = 

J Xjp 

2^ym^ -1- 2m^(2m^)^/^ 


3(^6^ 


xt^O 


3cy€A 


(26) 


where xt is the value of x at the top of the barrier, which is of the order of the size of the 
solitons, but we can approximate it as zero). This is the contribution to the action of the 
part of the instanton from the top of the barrier to the bounce point. 


C. Contribution of the interaction potential 

In this subsection we explain why the interaction potential does not contribute sig¬ 
nificantly to the energy and the action. If we imagine the sheep solitons separating, then 
between them the (j) field will be in its true vacuum but the "0 field will be in its false vacuum. 
Thus as they separate, the total energy will grow linearly. Thus the shepherd solitons give a 
compression that holds the sheep together. However as the shepherd solitons compress the 
sheep more and more, there will be an increase in the pressure, and the system will come to 
an equilibrium. There are many contributions to the pressure forces, the 0 kinks becoming 
steeper than their natural slope will add energy and hence pressure. But there will be a 
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contribution coming from the interaction term. Effectively the shepherd solitons would like 
to approach each other and annihilate completely in order to minimize their energy. The 
sheep solitons get in the way. The energy diminishes linearly as the separation of the shep¬ 
herd solitons goes to zero, with coefficient /3e^p. If one forces the shepherd solitons closer 
and closer, the value of the ip held at the point where 0 ~ 1/2, starts to move signihcantly 
away from ±a, the values for which the numerator of the interaction term vanishes. This 
will make the energy increase, to hrst approximation linearly, however with much greater 
coefficient > /3e^. Thus the dynamics can be simply modelled as given by a potential of 
the form 

Veff. = I3e^x - ^{x - xo), (27) 

where the second term is valid only for x ^ Xq and Xq is smaller than the size of the 
free solitons. For smaller values of x the repulsion becomes quite non-linear. However, 
the energy in the compression term is comparable to (Se^xo, which, and the corresponding 
pressure, for small can be made very small. This compression will be balanced by an 
equal, and hence small, but opposing, repulsive force. Thus this contribution to the energy 
of the solitons is arbitrarily small in the limit —)■ 0, which is what we assume. Therefore, 
the interaction term does not contribute signihcantly to the energy of the solitons. The 
most physical analogy can be made to nuclei, which have intrinsic mass energy, proportional 
to the number of nucleons, but which is modihed by the binding energy of the nuclei. For 
precision measurements the binding energies are important, but in the values of the overall 
masses, they are quite negligible. 


D. Contribution of the trajectory from the initial metastable state to the top of 
the barrier 

The evolution from the initial metastable conhguration to the top of the potential barrier, 
involves passage through conhgurations that cannot be well approximated by the free kinks. 
We would have to resort to numerical methods to calculate the contribution to the action of 
that part of the trajectory, however, we will in fact be able to conclude that it is negligible. 
With a simple parametrization of the conhguration in terms of hyperbolic tangents, we can 
get an estimate of the contribution to the action from the initial part of the trajectory. 
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FIG. 7. (Color online) (a) Metastable soliton configuration, (b) Solitons just touching, (c) Solitons 
separated by a macroscopic amount, (d) Solitons separated by a macroscopic amount. 

Taking 

il){x) = - (tanhf2(x + {^/2 + (2/a;))) — tanhI2(a; + {^/2 — (2/a;))) 

+ tanhr2(a: — (^/2 — 2/a;)) — tanhf2(x — (^/2 + 2/a;)) — 1) 

(j){x) = - (tanha;(x + ^/2) + tanha;(x — i/2 )), (28) 

this function is depicted in Fig. (7a, b). 

The numerically evaluated action for this trajectory gives an upper bound for the contri¬ 
bution to the action from the part of the trajectory which corresponds to the formation of 
the solitons at the top of the barrier from the initial metastable conhguration. The subse¬ 
quent contribution to the action has already been evaluated in the previous section in Eqn. 
(26). The numerical evaluation of the action is shown in Figure (8), with a zoom of the top 
of the barrier in Figure (9) in order to be conhdent that it is smooth. This behaviour of 
the potential is completely understandable in terms of the energy found in Eqn. (20), which 
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FIG. 8. (Color online) Potential for the parametrization of the bounce trajectory. 

corresponds to the part of the potential starting at the top of the barrier and descending 
with the very slightly sloped straight line. The bounce point is achieved aX xb = 2m^/ae0, 
which obviously behaves as o(l/e0). The trajectory from the metastable state to the top 
of the barrier, in the parametrized set of configurations, achieves the top of the barrier 
with a variation of x which is of the order of the size of the solitons, a size which is o(l) 
when compared with o(l/e<^). Using an analysis identical to that preceding Eqn. (26), the 
contribution to the Euclidean action, with a linear approximation to the potential, is given 

by 

r^T 

Se= dx^2m^x^m^ + 2m^ = + 2m^XT^^'^ ~ o(l). (29) 

Jo 

Thus it is perfectly reasonable to neglect this contribution to the Euclidean action. 
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FIG. 9. (Color online) Zoom of the top of the barrier of the potential for the parametrization of 
the bounce trajectory. 

IV. DISCUSSION 

We have computed the decay of a domain wall which traps the true vacuum inside the 
false vacuum. At hrst glance, one would think that such a domain wall would be classically 
unstable to dissociation, however it is easy to conceive of interactions which keep the soliton 
classically stable. The model presented here may be artificial; however, the general idea is 
clear. Any potential which prevents the passage of the sheep solitons through the shepherd 
solitons will do. 

There is furthermore no analogous surface energy that occurs in two [6, 7] and three [8] 
spatial dimensional examples, which could trap the true vacuum inside, at least classically. 
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However, in this paper we find that we can build on the analogy with the two spatial 
dimensional models [6, 7]. In those models, the true vacuum was separated from the false 
vacuum by a thin wall. The thin wall was achieved by having several topological quanta of 
one of the fields trapped inside the wall. The pressure of these quanta made the wall thin. 
It may well be that the thin wall monopole suggested in [8] in fact only occurs for large 
monopole charge as has been found in several large charged BPS monopoles [9]. 

In the present paper, we also trap N solitons of the field (j) inside the confining soliton of 
the field "0. The N solitons in principle exert a pressure and the overall size of the full soliton 
is proportional to N. The decay of the soliton can occur through tunnelling transitions. In 
the initial metastable configuration, the external false vacuum is due to the (j) field; the 
'll) field is in its true vacuum. However, inside, the V’ field goes to its false vacuum while 
the 0 field interpolates from the false vacuum to the true vacuum and back to the false 
vacuum multiple times, finally exiting the confining region in its false vacuum. Tunnelling 
transitions, essentially in the field 'ijj inside the confining soliton, will cut the soliton into two 
parts. The region between the two parts is in the true vacuum for both fields. This could 
happen at any point within the confining soliton where the 0 field crosses its true vacuum, 
which occurs multiple times if there are multiple 0 solitons confined. It will continue to 
occur until all the 0 solitons are free. The configuration will then resemble one half-soliton, 
linking the false vacuum on the outside to the true vacuum, followed by a string N of full 
solitons, in which the (j) field makes the transit from true to adjacent false to the subsequent 
true vacua while the 'ijj field makes the transition true to false and back to true vacua, 
followed by a final half-soliton interpolating from the true vacuum to the false vacuum on 
the outside at the other end. However, after any one transition, the false vacuum outside 
will be unstable to the groups of herded solitons from separating apart, converting false 
vacuum to true vacuum. The decay rate for each tunnelling transition is proportional to 




(30) 


where k is the usual determinant prefactor [10, 11]. We note that parametrically for a large 
enough, this decay rate can be unsuppressed. 
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